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ABSTRACT 

A simple classical mechanical model of the reactive scattering of a 
structureless atom A and a quasi-diatomic BC is developed which takes 
full advantage of energy, linear and angular momentum conservation relations 
but introduces a minimum of further assumptions. These are as follows: 

(1) the vibrational degree of freedom of the reactant (BC) and product 
(AB) molecules is suppressed, so the change in vibrational energy is simply a 
parameter; (2) straight-line trajectories are assumed outside of a reaction 
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"shell" of radius R ; (3) within this zone, momentum transfer occurs 

impulsively (essentially Instantaneously) following mass transfer; 

(4) the impulse I , which may be either positive or negative, is directed 

along the BC axis, which may, however, assume all orientations with 

respect to the incident relative velocity g . Assumption (3) is not as 

drastic as it appears since one may alternatively consider the reaction 

to have occurred on a potential surface such that the net force integrated 

over the entire "interaction region" yields the particular value of the 

impulse I . Differential and integral reaction cross sections are scaled 
2 

according to R (evaluated via molecular size considerations), but the 
model is not a member of the "hard-sphere" class. Potential surface in- 
formation enters the model in two ways; (1) by the choice of sign of the 
impulse, positive corresponding to net exit-channel repulsion ("late- 
downhill" Polanyi— type) , negative to a net approach-channel energy release 
("early-downhill"); (2) by the weighting function for impulse directions, 

I vs. g , l.e., the distribution function of (non-collinear) ABC config- 
urations. The model yields differential and total cross sections and 
product rotational energy distributions for a given collision exoergicity 
Q , or for any known distribution over Q . Numerical results are pre- 
sented for several prototype reactions whose dynamics have been well-studied 
in several molecular beam scattering laboratories. Many of the main experi- 
mental features are readily accounted for by the model without adjustment 
of parameters. However, it is difficult to reproduce the often--found 
"uncoupling" of the product angular and recoil energy distributions. 
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I. INTRODUCTION 

Molecular beam, chemiluminescence, and chemical laser experiments 

are rapidly expanding our experimental knowledge of reaction dynamics.^ 

Our theoretical understanding is also increasing, but it lags somewhat 

behind experiment. Although quantum, semiclassical and classical approaches 

2 

are all being intensively pursued, up to the present time classical tra- 

3 

jectory studies have probably contributed most to our understanding of 
the relationships between potential surface and the observable dynamics. 

The present paper considers a further-simplif ied, classical treatment of 
(model for) the 3-body exchange reaction. 

In light of advances at a more fundamental level one may question the 
role of such simplified models in the field of reaction dynamics. First, 
it is noted that when approximations are made to a more rigorous theory, 
in addition to facilitating calculation, one often tests the physical 
basis underlying the approximations and thus gains a deeper physical 
insight into the collision process. Second, one can hope to correlate 
trends for series of related reactions. Third, one can attempt to identify 
reaction features which can be attributed primarily to the existence of 
kinematic constraints (such as the conservation relations) so that those 
features most sensitive to potential surface characteristics may he isolated. 

In this paper a new classical model for reactive collisions of the 
type A + BC AB + C is developed, which is more general than the familiar 
hard-sphere models but which retains their attractive feature whereby 
all dynamical quantities follow from the equations of conservation of 
energy and momenta. The central theoretical construct of the model is 
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the impulse of a collision, defined as the force between two bodies integrated 
over the entire trajectory. This impulse I (units of momentum) may be 
viewed as the net momentum transfer in the course of the reaction. In 
some average sense, the impulse has extracted all the necessary information 
from the interaction potential. Because (or perhaps in spite) of this 
simplistic characterization of the interparticle interactions, the model 
applies both to reactions which occur on a single potential surface and 
to reactions involving potential surface crossings. 

Before proceeding, it is of interest to list the existing simple 
classical models. Excluded are statistical models,^ and those which deal 

5 

exclusively with ion-molecule reactions. 

Hirschfelder and Wigner^ and Hulbert and Hirschf elder ^ performed 

quantal calculations of transmission coefficients for particles under the 

influence of idealized potential surfaces, which were followed by classical 

8 

mechanical calculations for these same surfaces. The harpoon (electron- 

jump) model was developed by Magee^ to account for the large rates of 

flame reactions involving alkali atoms. One of the first fully classical 

models for reaction dynamics was the so-called billiard-ball (BB) model 

of Libby, applied to hot atom reactions to obtain the probability of 

reaction as a function of translational energy. BB models continued to 
XI 12 

be used until Cross and Wolfgang showed the inaccuracy of the BB model 
in predicting the energy-weighted total cross section (the so-called 
reaction integral). There have been criticisms^^ ’ of this conclusion 
and alternative BB models have been proposed. Most recently Baer and 
Amiel have extended the simple BB model to include idealized potential 
surfaces and have obtained improved results. 
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From the viewpoint of the classical kineticist the most important 

16 

model is that of Fowler and Guggenheim, i.e., the familiar hard-sphere 
model in which reaction occurs if the component of the relative velocity 
along the line-of-centers at collision exceeds a threshold value. Present 
placed this model on firmer theoretical grounds when he recast it in terms 
of the potential energy surface and its relation to activation energy. 
Extensions of this model have been made by Levine. Marron^^ has dis- 
cussed this model and a similar one which is based on consideration of 
the energy and momentum of the reactants and the products. Other modi- 
fications of the simple hard-sphere model for reactions include those by 

20 21 22 
Herschbach, Beuhler and Bernstein, and Grice. 

A different type of reaction model is the spectator stripping (SS) 

23 

one which was borrowed from the theory of direct nuclear reactions . 

This model was first applied to ion-molecule reactions and then to the 

2 ^ 20 

alkali atom-halogen reactions by Minturn ^ Herschbach has 

suggested a modification of the SS model which adds the internal momentum 

of BC to the final AB momentum to produce a broadening of the angular 

25 

distribution. Herman ^ have also suggested a modified form of the 

SS model, in which polarization forces are considered explicitly. Chang 
2 6 

and Light have extended it further by making additional dynamical 

27 28 

assumptions. Two other models, * very similar in nature, are directed 
primarily to the problem of estimating energy partitioning between trans- 
lational release and vibrational excitation of products. 

Another class of models, which leads to the present model, might be 

29a 

termed advanced hard-sphere models. Suplinskas has described a hard- 


17 
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sphere model for reactions A + BC, where each atom is represented by a 

29b 

sphere, and extended the model by including an attractive long-range 

force similar to that of Ref. 25. A final case is the DIPR model of 
30 

Kuntz et al. which has provided part of the Inspiration for this work,' 
The DIPR model was developed to fit the results of extensive 
classical trajectory calculations for reactions forming an ionic bond, 
and has been fairly successful when applied to ion-molecule reactions . 

The DIPR model allows calculation of the scattering angle by con- 
sidering the C atom which is deflected from its original path by a 
force along the original BC-axis direction. It may be viewed alterna- 
tively as a hard-sphere model in which A approaches BC with no inter- 
action until a fixed reaction distance is reached, then B changes 
partners and momentum transfer of a specified magnitude occurs instan- 
taneously in the BC-bond-axis direction. The rather drastic assumption 
of instantaneous momentum transfer may be softened by adopting a pseudo- 
hard-sphere view which defines the reaction radius as the point of mass 
transfer and defines the impulse (net momentum change) as occuring in 
some direction (e.g,, the original BC-axis direction) to be averaged over, 
and operating initially on the system A + BC and finally on AB + C . This 
is tantamount to enclosing the interaction region in a "black box" and 
spcifying the net momentum transfer. By observing how A and BC enter 
the box one can calculate the AB and C exit angles and recoil energy. 
If one assumes a particular separation at which mass transfer occurs, the 
rotational energy of AB may also be calculated. By varying the relative 
A— BC orientation and impact parameter, the angular and energy distribu- 
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tions may be computed. While Kuntz et al. have not described the DIPR 

32 

model in these same terms, the equivalence can be shown. A distinct 
advantage of the DIPR model over many others is that it is possible to 

30d ' 

derive analytical expressions for the energy and angular distributions. 

A second advantage of the DIPR model is that the free parameters ("loose 

screws") such as the magnitude of the impulse, orientation weighting, etc. 

30 33 

can be related roughly to potential parameters, ’ 

The present model, though closely related to the DIPR model, differs 

significantly in one way: instead of assuming the magnitude of the impulse. 

and calculating the final energy of the products, one specifies the latter 

and calculates the coresponding value of the impulse and thus the angular 

(and rotational state) distribution of the products. This is of some,. 

practical value in connection with modern molecular beam experiments 

34 

employing the velocity-selection and -analysis technique. The hope is 
that those features of the scattering which are due primarily to kinematic 
or conservation constraints can be isolated from those due to potential 
surface characteristics. When both the initial and final translational 
(c.m.) energies are specified the energy and momentum transfers are those 
required by the conservation equations. Calculations with the present 
model indicate that these constraints, together with a small amount of 
information about the potential surface, suffice to explain the main 
features of the angular distribution of products. 
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II. CONSERVATION EQUATIONS AND MODELS 

On the one hand a simple, easily-soluble model is desirable yet it 
must be sufficiently flexible to permit a description of a variety of 

26 

systems. Light's definition of a model is, "two trajectories connected 
by an assumption." Here the "trajectory" is that portion of a collision 
before (after) the reactants (products) feel significant mutual inter- 
action. The initial and final trajectories are then directly related by 
the conservation equations: 

+ c ~ ^ (1) 

/_ 4- j = U + (/ (2) 

-i- . (3) 

Equation (1) expresses conservation of linear momentum; A, B, and C 

represent the masses of the particles A, B and C and is the velocity 

of the i particle relative to the center of mass (c.m.). Equation 

(2) represents conservation of angular momentum; L is orbital angular 

momentum, and j is rotational angular momentum; primes denote final 

quantities. Equation (3) expresses conservation of energy; AD is 

o 

the appropriate difference in zero-point dissociation energies bond 
(positive for exoergic reactions) . In the present model there is no 
vibrational degree of freedom, so ^^ib^^vib simply a parameter. For 
purposes which will become clear, it is preferable to write the conservation 
equations differently, in terms of the impulse I ^nd the relative 
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velocity § = ~ ~C ‘ 

■1-^(6 + C)g^/(A+B+C) . Equation (1) becomes 


■A 'si-46 =- ^ + C w-c • 

Equating both sides to I gives: 

A£- g - CltlK- ?' = I (5) 

CA+^-^0 "" (A-b6b<^) ^ ^ ) 

or, alternatively 

j r = X + f ' «> 

where p = yg the initial relative momentum, p’ = y*g' the final 
momentum and f = C/(B + C). The partitioning of Eq. (1) to form Eq. (4) 
is somewhat arbitrary and could also be done: 

A 

which also yields Eq. - (6) with f = 1 . The difference is one of a scale 

30 

factor., Xhe first form has been chosen^ to agree with the DIPR model. 

Eq. (5) may also be written 


This implies that the particle C initially travelling with the velocity 
of BC , i.e., w is deflected from its path by the impulse I to 
velocity w . The scattering angle, say x > thus defined as the 
angle between p and p’ . The physical distinction between (4) and (4’) 
is that the impulse occurs after mass transfer for (4) but before mass 
transfer for (4’). 
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The conservation equations are not, in themselves, sufficient to 
determine the collision dynamics, i.e., given a set of initial conditions 
the final trajectory cannot be determined. This is due to the fact that 
there are 7 equations, (two vector, Eqs. (2), (6) and one scalar, Eq. (3.)) 
and 10 unknowns: L', j', g’, and • The model is a set of 

assumptions which implies values for any three of these ten unknowns. 
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III. DERIVATION OF THE MODEL 

Initial experimentation with various hard-sphere models revealed 

that the requirement of momentum transfer along the line-of-centers 

35 

produces only "backscattering" of the reaction products. Figure la 
illustrates the vector relationship expressed in Eq. (6) for conservation 
of momentum; it can be seen that if the impulse lies in the same (forward) 
direction as the initial relative velocity, as it almost always does for 
line-of-centers momentum transfer, the final relative velocity will be 
directed backward. In order to obtain a model that will produce forward 
and/or backward scattering the impulse must be permitted to assume all 
angles with respect to the incoming p . One method of achieving this 
goal is to discard the line-of-centers assumption for I and assume that 
I must be directed along the B-C bond axis. Since the orientation of the 
target molecule BC is not restricted, I can then take on all directions. 
The concept of a hard spheres is used only to determine the "reaction 
radius" R (the distance between A and the c.m. of BC) which is the 
point on the initial trajectory at which mass transfer occurs. It is 
assumed that momentum transfer occurs impulsively, though not as it would 
for hard spheres. The assumption of instantaneous momentum transfer serves 
to simplify the solution of the equations of motion in that the conserva- 
tion equations then contain all the system dynamics. If one prefers, the 
reaction may be viewed as having occured on some potential surface with a 
net force, integrated over the entire interaction region, in the direction 

of I and having a magnitude |l] . These two viewpoints are formally 

36 


equivalent. 
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The coordinate system is shown in Fig. lb. The direction of the 
z-axis is that of the initial relative velocity g ; the orientation of 
the target upon collision is given by the angles 6 (the rotation of C 
about the BC center-of-mass c.m. in the XZ plane) and (j) , the azimuthal 
rotation about the z-axis. The scattering angle x is defined by 

X f f fO » 

Squaring Eq. (6) gives 

4 j f f . (8) 

Defining the dimensionless parameter 

(X ^ f V ff ^ (9) 

Eq, (8) is rewritten 

- C*" + ^ 

The Z“Component of Eq. (6) is 

= i ■ (U) 

Then 0 can be expressed in terms of x • 

= ( I odCAh<5x)/( l --Zo^OK>y (12) 

Eq. (11) may be substituted into (10) to yield the inverse relationship 
for x(0) : 


( 12 ') 
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where the negative square root is implied if a ^ 1 ; both roots are taken 
if a < 1 and 0 _< sin0 _< a • Note that the scattering angle x is 
independent of the azimuthal angle cj) of the target. It can also be shown 
that the azimuthal angle of the trajectory, say ip , is identical to ^ . 

Eq. (12') indicates that the scattering angle x is independent of impact 
parameter, b (the perpendicular distance between A and the z-axis at 
collision). The final rotational energy, however, depends on both 

b and \p . 

A. Potential Surface Information 

Before going on to derive expressions for product angular and rota- 
tional energy distributions it is of interest to note that qualitative 

potential surface information may be introduced at this point by two means. 

First, the sign of the impulse can be chosen as positive or negative. 

37 

Following Polanyi, this might correspond to whether the potential is of 

the ” early-downhill" or "late-downhill" type, i.e., whether energy is 

released primarily in the approach channel, A + BC, or the exit channel, 

AB + C. Since in the model the impulse occurs after mass transfer (i.e., 

Eq, (4) instead of (4')), a positive impulse corresponds to a net exit- 

channel repulsion while a negative impulse indicates a net approach-channel 
38 

energy release. These concepts are discussed ia more detail in Ref. 30. 

A second means of introducing potential information is to weight 

different impulse directions, i.e., target orientations, differently. 

21 

For example it is known that in the reaction M + CH«I MI + CH., the 
alkali atom attacks the iodine end of CH^I preferentially. On the other hand. 
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/for an alkali 'attack on HBr . (where the Br- atom is much larger than the 
H atom) much less :orientation dependence is .expected... .The simplest 
possible weighting, function is one /involving only the. angle* 0 .(see 
Fig. lb).> A convenient functional form .is., . . .... > 




(13) 


where, A .is a normalization constant. The plus or minus, sign, is chosen 
according to whether the B or the C end of the molecule is reactive. 

The power n determines how rapidly drops to zero. Note also that 

by virtue of Eq. (11) and because of the symmetry in the form of Rg ^ 


changing the sign in front of cos0 is equivalent to changing the sign 
of the impulse I . A sign change may thus be considered either a change 
in which atom of BC is most reactive o^ a change in the sense of the 
impulse I . (A change in both leaves the sign in (13) unchanged.) This 
equivalence exists only because of the symmetry inherent in Rg due to 
the fact that is assumed to be directed along the BC axis. 

The weighting function (13) is not. simply related to the A-B-C 
configuration angle, say C > i-e., the angle between R and the BC 


axis (see Fig. lb) 


. Another choice 


, f or^ the,, weighting . ^func tiqu..^i^^ 


’>:l J 


where. A’ 

A A nC 'i ' L. 




i£ i ' < i ‘ ' ' A' ’ ;f V. 

is a constant and cosC 


is given by 


U4) 






T S-Av © ’ 4- c os © ' £<>-0 Y I 




'(15) 
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The angle ? depends on both the orientation (0,(f)) of BC at collision 
and the impact parameter b = R sin Y , where Y Is the angle between 
R and g (see Fig. lb). Averaging over impact parameter must now be 
carried out explicitly. For the case n = 1, averaging over cj) and b 
yields ; 




SwuY CiX K 


1 ± L C4>d e 


(16) 


(Similar equations can easily be found for n ^ 1.) Comparing Rq and 

<T„> we see that the effect of weighting by cos^ instead of cos0 

y 

is to give a higher relative weight to target orientations in the region of 9 

near tt/ 2. This is reasonable because both and give the 

u y 

largest weighting to collinear configuration but fails to take 

y 

account of the fact that collinear configurations exist for all 6 

In both of the two weighting schemes above a step function ^ (R-b) 

has been assumed for the impact parameter dependence is the Heavi- 

39 

side function). A more realistic form for P, (b) has a maximum at 

D 

b = 0 and decreases monotonically to zero at some b , say R . This 
form reflects the fact that collisions at large b have only a small 
radial velocity component; even a collinear collision at large b is 
less likely to result in reaction than a collinear collision at b = 0 . 

As an illustrative calculation the functional form P^acos y has been 
used. The major effect of including this P.^ is to counteract the change 
noted in going from Rg to <T^> , that is target orientations with 
0 near 0 and ir are reemphasized over orientations near 6 = tt/2 . 

The weighting functions Rg , <Tg> and ^Tg*P^> are plotted in Fig. 2 
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for the chosen P^«:qos y and n = 1 . The brackets indicate averaging 

over b and (j) . Because of the similarities between Rg and the most 

"realistic" weighting function <T^*P,> and because there is no simple 

y D 

a priori way to determine either T^ or P, , the convenient, simplest 

U b 

form, Rq , is employed for the remainder of the discussion. 


B. Differential Reaction Cross Sections 

The formula for the differential reaction cross section may be 

40 

derived from a formal definition of the cross section in terms of 
the model variables 






(17) 


where T] refers to the pair of variables b,3 (3 is the azimuthal 

angle of R ) , T refers to the orientation 0,c() (see Fig. 1) and 
0 ) to the scattering angles • The differential "area" on the 

left is proportional to the product of probabilities that for collisions 
in the range ri to n + dp and for target (BC) orientations in T to 
T + dT , the AB product will be scattered into a solid angle in the 
range w to o) + d^^ . The details of the derivation may be found in 

Appendix I. For the present model, Eq. (17) reduces to 




oc" 


(^oc-e 







(18) 


where O = ttR and Pg is the desired configuration weighting function. 
The argument is subscripted to indicate the one-to-one relationship be- 
tween X 0 . Eq. (18) will be discussed further in Sec. IV and 

applied to a number of systems of experimental interest in Sec. V.l. 
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C. Total Reaction Cross Section 


The reaction cross section O is found by integrating Eq. (18) 

K 

over all solid angles. Because of the absolute value in the numerator 
of Eq, (18) two cases must be distinguished: a > 1, a 1 . Assuming 

that is the step function (R-b) and Pg = 1 , direct integra- 

tion yields 


Cr r: 


CoC > 1 ) 


If a weighting function other than (R-b) ^ is used for impact parameter 

and/or if any weighting function is introduced for orientation, the total 

reaction cross section is necessarily diminished. The actual decrease 

depends on the explicit form chosen for P^ and/or P- . Equation (19) 

b b 

19 

is in agreement with the previous work on simple hard spheres where it 
was found that the total reaction cross section has an energy dependence 
only when a ^ 1 . Note that the present a differs by a factor of 
C/ (B + C) from that in Ref. 19. 

D. Product Rotational Energy 


The rotational energy, products may be determined 

via Eq. (2); details are given in Appendix II. It is found that 


depends 




( 20 ) 
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upon the system parameters R, y and y' , g and g' (which are 
given, the direction g' defined by the scattering angle), and b , 0 , 

(f) defining the reaction configuration. One must average over b , (p 
and ip (including any desired weighting function) to determine the con- 
tribution from all configurations, to yield the quantity 
brackets < > denote such an average over b , cj) , and \jj . This 

calculation has been carried out in Appendix II. Then “^^rot^ depends 
only on y the ratio g'/g in addition to the system parameters. 

The result of this averaging can be expressed by means of the dimension- 
less quantity 6 , defined 

^ rot 



where f = C(B+C) Equation (21) is applied to several systems in 

Sec. V.2. 
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IV. BEHAVIOR OF DIFFERENTIAL REACTION CROSS SECTION 

The model predictions depend on the choice of values for R , for 
P(a) (or equivalently P(E’)), i.e., the distribution of a ; the sign 
of the impulse; and the configuration weighting function. In this section 
the behavior of the differential cross section is examined as a function 
of each of these quantities. 

The reaction distance R enters the differential cross section for 

the reaction first as a scale factor, i.e., scaling the differential and 

2 

total cross section according to R , and secondly as it affects the 
final rotational energy distribution (cf. Appendix II). Insofar as 
angular distributions are concerned the effect of R may be ignored. 

The value of a and its probability density depend upon the final 
translational energy and its distribution, P(E’) . In the present model 
these quantities are either to be guessed or taken from experiment. It 
will be seen that the angular distribution calculated from the model is 
not extremely sensitive to the choice of P(E*) . Unfortunately, this 
implies that the inverse use of the model to determine P(E’) from the 
angular distribution is unlikely to yield useful results. 

The final two quantities, the sign of the impulse and the weighting 
function, serve to characterize the potential surface for the reaction. 
Here one must exercise chemical intuition and allow the experiments to 
serve as confirmation. 

2 2 

In Fig. 3 angular distributions d CT (x)/d co are plotted for various 

K 

choices of alpha and the weighting function Rq given in Eq, (13). As 
mentioned before, because of the inherent symmetry of R^ 


a change of 
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sign is equivalent to either interchanging the "reactive" end of BC or 
changing the sign of the impulse. If the limit of Eq. (18) is taken 
as a we obtain 


(JO hit 

Thus, for large a and no weighting (n = 0) the angular distribution 
is nearly isotropic. When weighting is included, i.e., ^ 1 , the 

o 

angular distribution reflects entirely the weighting function. On the 
basis of this model, we expect that for experimental situations char- 
acterized by large a , the angular distribution of the products will 


41 

be very sensitive to potential surface features . This condition obtains 

in either of two situations: first, when the collisional exothermicity 

Q = E’ - E is very large (i.e., the AB molecules are produced in low- 

lying internal states and AD^ is large), and second when the mass B 

is much larger than the mass of A or C . This last result is easily 

B B 

seen by considering the limit of a for y , -^ > > 1 . Here 




r ^ \ y. 


As a approaches 1 from above, the unweighted (n = 0) distribution 
increases rapidly in the forward direction. When a weighting function 
is included, the angular distribution can be shifted either in the for- 
ward or backward direction or even ( e. g . , in the case a = 2 and n = + 1) 
to produce side-scattering. Forward scattering obtains when a 1 
since a large fraction of target configurations yield products in the 
hemisphere; in fact, all configurations with 0 ^ 90® yield products 
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within a C£. 1° forward cone. Inclusion of a weighting function serves 
either to emphasize the minority configurations which produce proportion- 
ally more products in the backward direction or vice versa. If the 
magnitude of the impulse is computed for a value of a just slightly 
greater than 1.0, say a = 1.001, it is found to vary between 2.0 and 
0.0 (in units of the initial momentum yg ) with I ^ 0.01 for those 
configurations which yield foward scattering. 

It is instructive to compare the a = 1 case with the predictions 

20 24 

of the spectator stripping (SS) model. * The essence of the SS model 
is that atom C is a spectator to the reaction and will continue with 
velocity and direction unchanged after loss of B . Thus Wg^ = w^ , 
which implies together with Eqs. (4) and (5) that 1=0 or that 

h = f' . . 

Using Eq. (7) the scattering angle may be computed, which yields the well- 
known result that the products are always scattered at X “ 0 • The SS 
limit cannot be obtained directly from Eq. (18) because we can only 
include one of the three SS constraints exhibited in Eq. (24), that is 
a = 1 . The limit (a ->■ 1) applied to Eq. (18) yields an infinite differ- 
ential cross section for x 0“ > reminiscent of the classical divergence 
in the differential elastic scattering cross section at 0° for all but 
rigid spheres. The SS limit can of course be derived directly from 
Eq. (6) where all three constraints may be included. 

For values of a < 1, Eq. (18) has a zero at x “ cos a . The 
physical reason for the zero is that not all possible reaction configura- 
tions can produce a reaction. Fig. 4 is a plot, for a = 0.8, of the configur^' 
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tion angle 9 versus the scattering angle x • configurations 

0 > 53°, reaction cannot occur satisfying the conservation relations, 

This condition is not unique to this model but is found in other models 
in which energy partitioning is specified. When the angular distribu- 
tion is averaged over any reasonable final energy distribution the zero 
disappears and predominantly forward scattering obtains. Inclusions of 
the weighting function, Rg , will modify the angular distribution; but 
forward scattering persists. As a tends to zero the differential cross 
section and the total cross section (cf. Eq. (19)) also tend to zero. 

The situation characterized by small a arises only for instances of 
negative Q . 

It should be noted that the joint angular -recoil energy distribution 
is a coupled function of x ^nd E' (via a ) . For "rebound" scat- 
tering (i.e., positive impulse) the angular distribution shifts from 
sharply backward to broad, as alpha decreases. When the impulse is 
negative, the angular distribution is found to be broad, slightly favoring 
the forward direction for large a , and as a is decreased the scat- 
tering is progressively concentrated in the forward hemisphere. 



21 


V. APPLICATION OF MODEL 
A, Angular Distributions 

In this section the differential reactive cross sections are cal- 
culated for several systems which have been studied experimentally. It 
will be seen that for reasonable choices of the weighting function, good 
agreement can be obtained between the model and experimental results. 
Further, the model can explain the observed trends of the reactive scat- 
tering for members of a homologous series. The results described show 
thereby that many of the gross features in the differential cross section 
may be attributed entirely to constraints of conservation of energy and 
momenta. A posteriori at least, the feasibility of estimating appropriate 
weighting functions is demonstrated. 

Rb + CH^I 

The reaction of an alkali with CH^I is an example of a repulsive 
interaction exhibiting backward ("rebound") scattering in which a sub^ 
stantial fraction of the reaction exothermicity is converted into 
translational energy. 

21 

Beuhler and Bernstein have studied the reactive S3rmmetry for this 
system caused by orienting the methyl iodide molecules. They find the 
Rb-ICH^ orientation to be at least five times more reactive than the 
Rb-H^CI orientation. While we cannot make explicit use of this data 
to determine a weighting function, these results suggest a choice of n 
in the range 1-2 for the weighting function Rq . Recent experiments^^ 
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indicate that the most probable E' is large, probably about 15 keal/ 

aXXS 

mole, which for Rb + ^^ 3 ^ gives a most probable a of 15 . Fig. 5 gives 

the results for n = 1 , positive impulse (repulsive interaction), for a 

in the range 10 to 20. The calculated curves resemble the published 

45 

e?cperimental results for this system. Because of the high value of 

a , the angular distributions are nearly independent of a . This result 

44 

is in accord with the experimental finding that the angular distribu- 

46 

tipn is only weakly coupled to the energy distribution. 

M + X 2 “ Trends 


In this section the trends expected in the alkali halide systems 

are examined in order to compare with those observed in the angular 

distributions for K + CI 2 » Br 2 and I 2 by Grice and Empedocles^^ and for 
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Li, Na, K, Rb, Cs + Br 2 by Parrish and Herm. 

For the X 2 systems we may use the model in its simplest form ( i. e. , 
employing a minimum amount of experimental information) . Three parameters 
must be assigned: the sign of the impulse, the value of a and index 

n of the weighting function Rg . The alkali-halogen reactions are 
known to be of the ’’attractive** type as evidenced by their extremely 
large cross sections, thus implying a negative impulse. The simplest 
assumption regarding the value of the most probable a is that the 
reactions are ’’collisionally thermoneutral," i.e. , Q = 0, so that a 
is determined by mass factors alone. The resulting values of a(Q ^ 0) 
for K + Cl 2 > Br 2 and I 2 may be found in Table 2. 
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The choice of weighting is only slightly more difficult. Some weighting 
is clearly indicated since configurations such as A-C-B, where A departs 
with B, should have considerably less weight than A-B-C. (Recall that 
even though B and C are here identical X atoms, the model distinguishes 
between B and C.) The form of the weighting function Rq , i.e., the 
choice of the weighting parameter n is however, not extremely crucial, 
and a choice of n = 1 seems reasonable. 

With these values for the model parameters the calculated angular 
distributions are shown in Fig. 6. They agree well with the experimental 
results of Ref. 47. 

Since more detailed data exist for K + Br 2 and K + it is worth- 
while to apply the model in its more flexible form to determine the 
"experimental" values for the sign of the impulse and weighting parameter 

by fitting the observed angular distributions using experimental values 

for a . For K + Br^ the attempt is made to fit the angular distribution 
49 

of Birely et al . , folding in the product recoil energy distribution 
function of Warnock et al .^^ The product energy distribution is con- 
verted to a distribution over a at the nominal collision energy E . 

The optimum, fitted, value for the weighting parameter was found to be 
n = 0.7, with a negative impulse as anticipated. The differential cross 
section averaged over a is not extremely sensitive to n , but is, of 
course, very sensitive to the sign of the impulse. (Switching from nega- 

51 

tive to positive impulse shifts the products into the backward hemisphere. ) 


For K -h complete flux-velocity -angle contour maps 

52 

as a function of angle and velocity) are available 


(product intensity 
for several E 

tr 
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The optimum parameters found by fitting the entire flux maps turn out 
to be the same as the above values for the K + Br^ system. 

A negative impulse and a weighting parameter n = 0.7 are used for 

the Li, Na, K, Rb, Cs + Br^ systems. The values of a used are those 

listed in Table 2. These calculations give similar angular distributions 

to those of Fig. 6. Table 1 lists the results as X 2/2 ’ width of 

2 2 

the angular distributions d CT^(x)/d co at half-height, compared with 
the experimental ^j^/2 The trend observed by Parrish and Herm^^ 

is found, supporting their assertion that the broadened LiBr angular 
distribution may be attributed to a "mass effect in the overall reaction 
dynamics . " 

Note that the present model has nothing to say about how the avail- 
able energy is partitioned between product excitation and recoil energy, 
but givsn this partitioning, the kinematics require a broader angular 
distribution because of the small Li mass which results in a large value 
for a . Table 1 also shows that as the alkali mass increases the angular 
distribution becomes more peaked in the forward direction. However, the 
observed K, Rb, Cs + Br 2 angular distributions^^ appear to be virtually 
identical. Possibly the weighting index n = 0.7 should be decreased 
slightly for the larger atoms to permit more reaction configurations. 

This change would bring the calculated Xj^/2 for the K, Rb, 

and Cs systems closer together and have little effect on the Na and Li 


systems. 
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Table 1. Width at half-height, X ]_/2 » angular distributions, 

2 2 

d o_(x)/d to , for the M + Br„ reactions (peaked at X ~ 0) 
R ^ 


M 

Xi/2 (deg.) 

Model^ 

Calc. 

Exptl. 

Li 

85 

100 

Na 

70 

70 

K 

60 

60 

Rb 

55 

60 

Cs 

35 

60 


Model parameters: negative impulse, weighting index n = 0-7, 

a values from Table 2 . 

K + HBr 

54 

The systems K + HBr, DBr were studied by Gillen et al . Because 
of the pathological kinematics for these systems the derived c.m. 
angular distributions are quite uncertain and a meaningful determination 
of the model parameters by a fitting procedure cannot be made. One can 
say, however, that because the alphas are large [at? 11(K, HBr), 

8 (K, DBr), and 7 (K, TBr) ] , very broad distributions are predicted 
unless the weighting parameter n is much greater than unity. (This 
seems unlikely considering the size of K and Br compared to hydrogen.) 
Angular distributions obtained from the model for any reasonable choice 
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of weights (with either positive or negative impulse) are found to agree 
qualitatively with the HBr and DBr data of Ref. 54, but predict that the 
KBr from the reaction with TBr should be peaked slightly more forward 
than for K + HBr or K + DBr. Unfortunately, this is contrary to the 
findings of Martin and Kinsey. From Ref. 54 the energy distributions 
for the H- and K- systems are both broad, with the D- distribution 
slightly broader than that for H . If there were a considerable contri- 
bution to the T- distribution from higher E' ( i.e . , higher a), this 
discrepancy could be explained. 


B . Product Rotational Energy Distributions 


Based on the analysis given in Appendix II, calculations have been 

made of the angular distribution of the average change in rotational 

angular momentum of product AB with respect to reactant BC. The quantity 

given in Eq. (21) has been calculated for various scattering angles 

X . For the systems chosen the values of the parameters r , R (and 

BC 

g'/g) used are those listed in Table 3. 


For the case of K + HBr, TBr, the calculations indicated that for 
both systems only weakly dependent on angle, declining with 

X smoothly from a value of 1.1 at X = 0° to 0.9 at x = 180°. Averaging 
over angles, assuming an isotropic angular distribution of products, 
yields a quantity which is related to tbe average finql rota- 
tional energy by the following equation: 


E 


roi 




< s > 


^7T 


+ £- 


ret 


' ( 25 ) 
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Here 


and are, respectively, moments of inertia of BC and 


AB . For E ^ 0.15 kcal/mole, E’ turns out to be 8.4 kcal/mole 

rot ^ rot 

for K + HBr and 1.8 kcal/mole for K + TBr . For K 4- HBr, 8.4 kcal/mole 

57 

exceeds the energy available for rotational excitation. Equation (25) 

may be used to compute product rotational energy at specified angles by 

substituting 6 for <6 . Examination of the K + HBr results 

at each angle reveals that exceeds 6.2 kcal/mole at every angle. 

The difficulty can be traced to the extremely sensitive dependence of 

E* on the value of R while 6 is only very weakly dependent on 

rot rot 

the choice of R . A change in R from 3.34 to 2.5 modifies the 

vs. X plot imperceptibly but will reduce E’ to less than 6.2 kcal/ 

rot 

mole at all angles “ 1.40 corresponds to “ 6 .2 kcal/mole) . 

Note, however, that a large value of E’ for K + HBr is expected be- 

rot 
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cause the mass ratio demands that L = j’ . 

- Because of this extremely critical dependence on R , which arises 
in converting angular momentum to rotational energy or vice versa, one 
attractive potential application of the model must be set aside. While 
general features of the distribution (flat, peaked forward, etc.) 

are significant, no real significance may be attached to the absolute 
values of • This precludes the possibility of invoking the con- 

servation of energy to obtain constraints on the angular distribution of 
products, based on rotational excitation of the products. 

Figure 7 presents the results X systems Li, 

K, Cs + Br^ . The behavior seen there is in general agreement with 
30 

Kuntz et al . who found a positive correlation between scattering angle 
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and enhancement of product rotational energy in their potential surface 

calculations for K + Br 2 • These authors do not report a plot such as 

Fig. 7, but only a single number which is the correlation coefficient 

for rotational excitation and scattering angle. Our results indicate 

that the excitation function vs. x broad, slightly favoring the 

forward direction. This prediction is actually in closer agreement with 
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the experimental study by Grice et. al. , who found that the reactively 

scattered RbBr from Rb + Br 2 could be characterized by the same rotational 

temperature at all angles for which they made observations. The most 

outstanding feature of Fig. 7 is the bimodal excitation function for 

LiBr. The bimodal distribution persists for a wide range of (gVg) 

values and for 0.2 < r /R < 1. (Values of r_^/R > 1 were not investi- 

rJL — BC 

gated.) > computed from Eq. (25), for Li, K, Cs + Br 2 is 0.5, 1.2, 

and 2.5 kcal/mole, respectively, with = 0.6 kcal/mole. If the 

angular distributions are taken into account in the average and consider- 
ing that for M + (B) reported in Table 1 of Ref. 48 are probably 

too low, good agreement between the model and experiment^^ ’ can be 


obtained . 
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Table 2* Summary of relevant experimental results from the literature. 



_>• 

E(kcal/mol) 

E’ (kcal/mol) 
mp 

u/y' 

a 


Rb + CH^I ^ 

1.5 

~ 15 

3.81 

15.3 

90 

Li + ^ 

2.0 

6.6 

0.160 

9.09 

100 

Na + Br ^ 
2 

1.28 

5.2 

0.447 

6.03 

70 

K + Br^ ® 

1.23 

3.7 

0.657 

4.28 (2.47^) 

60 

Rb + Br ® 
2 

1.01 

2.8 

1.03 

3.28 

60 

Cs + Br^ ® 

1.06 

1.8 

1.25 

2.33 

60 

K + Cl S 
2 

— 

— 

1.05 

1.95 ^ 

— 

K+ I 2 « 

— 

— 

0.471 

f 

2.91 

— 

K + HBr 

2.8 

1.5 

26.6 

11.5 


K + DBr ® 

2.8 ^ 

1.5 ^ 

13.5 

8.15 


K + TBr ® 

2.8 ^ 

1.5 ^ 

9.08 

6.70 



a. Nominal value of average initial relative energy. 

b. Estimated from unpublished results on K -f CH^I by A. M. Rulis and 
R. B. Bernstein 

c. Ref. 48. 

d. Ref. 61. 

e. Ref. 49. 

f. Alpha computed assiopiing Q = 0 . 

g. Data not entered if not referred to in text. 

h. Ref. 54. 


i. Assumed identical to K + HBr 
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Table 3. Summary of parameters used in the calculations. 




(g'/g) ^ 

mp 


02 

CrCa^) 

o 

R(A) 

K + 

HBr 

3.77 

1.41 

35 ^ 

3.34 

K + 

TBr 

2.21 

1.41 

35 ^ 

3.34 

U + Br^ 

0.726 

2.3 

115-146 ® 

~ 6.4 

K + 

Br^ 

1.41 

2.3 

220-260 ^ 

~ 8.7 

Cs + Br 2 

1.46 

2.3 

370-380 ^ 

~ 11 

a. 

Computed by (g*/g) 

mp 


The error 

resulting from 


using 

E instead of 

mp 

g is probably less than 

mp 

the uncertainties 


in the 

energies . 




b. 

Ref. 62 

• 




c. 

Ref. 54 

• 




d. 

Assumed 

equal to that 

for K + HBr. This is in disagreement ylth 


expectations based on 

the K + HBr, DBr 

experiments 

where a signi- 


ficant 

difference in 

was found. 

However, it 

is believed that 


the difference observed is largely due to the quantal Isotope effects 
which cannot properly be treated by the present model. 


e. Ref. 48. 


f. Ref. 49. 
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VI. CONCLUSIONS 

A simple kinematic model for reactive collisions A + BC AB + C 
has been developed which includes aspects of mass and momentum transfer 
and for which the equations of motion are readily soluble. Required as 
input are the masses of the particles, the initial and final transla- 
tional energies, and a value for R , the distance between A and the c.,m. 
of BC at which mass transfer occurs. Information on the potential 
surface is introduced via a weighting function for the direction of net 
momentum transfer, or alternatively reaction configurations, and by the 
sign (positive or negative) of the net momentum transfer. The model is 
formally equivalent to an impulsive model, yet because of the definition 
of net momentum transfer or impulse, it has greater generality. Attractive 
and repulsive systems may be treated. The model admits an analytical 
solution for the differential and total cross sections. No internal 
degrees of freedom are included so vibrational energy enters only as a 
parameter through the conservation of energy equation. 

The differential reaction cross section exhibits a wide variation 
in shape depending on parameter values. By changing the initial and 
final energies, the masses, and the configuration weighting, it is pos- 
sible to obtain broad or narrow distributions peaked forward, backwards 
(or even sidewards). The reactive scattering for Rb + CH^I; M + X2 ; 
and K + isotopic HBr molecules has been examined to determine if reasomble 
choices for the potential parameters can be made based on gross, "chemical" 


considerations . 



32 


One purpose in developing this model was to isolate those features 
of the differential reaction cross section which could be attributed 
mainly to kinematic constraints. To do this with a single model, a 
certain amount of adaptability must be added to the model by including 
parameters characteristic of the potential. The Rb + example dem- 

onstrates that the model is applicable to repulsive as well as attractive 
systems. The model predicts only very weak coupling between the recoil 
energy angular distributions, at least when a is large 

(as it is for Rb + ^ series) . For a near unity, 

significant coupling is predicted. For the M + X2 series the trend in the 
angular distributions for K + Cl2> Br^ > ^.nd has been established using 
the model in its crudest form. The assignment of the weighting parameter 
n was refined by choosing one system in the series (K + Br^) as the 
touchstone and fitting the experimental results to optimize the parameter, 
and then predicting fairly well the trend found experimentally in the 

Li, Na, K, Rb, and Cs + Br^ series. These calculations also confirm 

48 

Parrish and Herm’s "mass" effect for Li. For the K + HBr series it 
is found, with little or no potential input, that a broad angular dis- 
tribution of products is to be expected on kinematic grounds alone. 

Generalized kinematic models should prove valuable in two respects: 
first, by providing a means for isolating features of reactive scattering 
which are dictated primarily by kinematics and thus allowing one to con- 
centrate on those features directly related to particle interactions, and 
second, by providing a basis for understanding or predicting trends in a 
series of reactions involving similar species. Another possible applica- 
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tion would be to provide a convenient source for rapidlyr-calculated 
reactive trajectories which might be used in a s^iclassical treatment 
of reactive scattering. 
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appendix I: DERIVATION OF DIFFERENTIAL REACTION CROSS SECTION 

In general terms one may view the "generalized differential reaction 
cross section" as a probability density f (x) which is a function of a 
multidimensional variable x , and which is related to a distribution 
function F(x) by f(x) = dF(x)/dx . The quantity dF then is the 
probability that the variable assumes a value in the range x and x + dx* 
For the reactive scattering problem x subsumes the variables b and 
3 which are the impact parameter and azimuthal angle; 9 and (\> which 
are the target orientation angles; and y ^ which are the scat- 

tering angles. The infinitesimal element of probability is a six^ 
dimensional quantity 

aVOt.T,^) = W a^(iC-c)ci"H 

/[ ^ 

CO 

where S, G, and H are two-dimensional distribution functions for the 
appropriate variables. In writing this equation it is implicitly assumed 
that the variables ri > T , and cj are independent, which is not 
the case in general. However, quite often the coupling between sets is 
simple, so one can proceed with the independent formulation for clarity, 
incorporating coupling as it arises. It is worth emphasizing that all 
the dynamical information contained in the model pertinent to the scat- 
tering angle must be included in the distribution functions. 
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As an example consider the simple hard sphere (SHS) model in which 
T is a null set and b and X related uniquely. Equation (I.l) 

is simplified: 


Since b and 3 are to be integrated, we consider x ^ depen- 

dent variables. The function S is not a distribution function in the 

2 

usual sense since it is "extensive" with a maximum value of 27TR which 
2 

comes from b 3 . (This distribution is therefore not normalized to one 

The differential cross section for this model is 



The P functions are weighting functions which contain the rest of the 

model information not explicitly contained in (1.3). In this case they 

are all unity except for , which is a step function with value 1 for 

b £ R and zero for b > R . The delta functions arise because the 

relationship X]^ “ 2 cos ^(b/R) establishes a one-to-one correspondence 

between x ^ ij; = 3 because all forces are impulsive. The 

primed b’s are the problem variables, the unprimed b’s are functional 

(dummy) variables. The normalization for H is determined by 

N//6(iJ;- 3)6 (x^ “ X|^»)sinx^dxdi|; = 1 , where b’ and 3 are independent 

variables, distribution function H may also be written such that 

2 

d H(o)) = 6(i|j-3)5 [cosx^ - cosx^? ]d(cosx^)di|; . 
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The integration over 3 is easily performed and serves to remove 
the function 6 . To perform the integration over b’ the inde- 
pendent variable must be made explicit in the function • 

This can be accomplished by using the identity 

^IgCx)] = ^ Is' (x^) » where are the roots g(x^) = 0 . 

For this model there is a single root b = b' and Eq. (1.3) becomes 


which gives the well-known result 

cf^(T C^) - ^ 

or 

— _b I Jl b I ^ 

cft2-U0 ( d% I 

The subscript b has been suppressed since it is now redundant. Equation 

(1.6) is the well-known SHS result. Upon inserting the relation be- 

2 2 2 

tween b and x > one finds d O /d oo = R /4 . 

brlb 

The general result, including G(T) , is 


I SUa^ 

cix I 


cUif 


(1.5) 


( 1^)1 



(1.7) 


where the P’s must be determined from the model under study. The order 
of Eq. (1.7) may be increased by including distribution functions for 
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translational and/or internal energy. In keeping with the goal of de- 
riving a simple model these possibilities are ignored. One must examine 
next the coupling between rj , T and go . In the present model, the 
probability that a certain reaction occurs yielding products in go to 
CO + dGO when the target has an orientation in T to T 4- dT ^ does not 
depend on impact parameter. Separability does not obtain when configura- 
tion weighting functions are used which depend on the impact parameter. 
This situation is discussed in the section on configuration weighting 
where the net result of including such weighting is to produce a constant 
factor which may be included in G(T) . Coupling between 00 and T 
will always exist for any model in which target orientation plays a role 
in the angular distribution of products. In the present deterministic 
model there exists a one-to-one relation between (0,c()) and (X>4^) which 
may be expressed as delta functions ^(Xq^Xqi) 6(cj)-i|j) . The theta 

subscript on X0 indicates the choice of y dependent variable 

where cosXg is given by Eq. (12^ in the text. Equation (1.7) may 
be expanded as follows 



where 


P,, can be taken as 1, may be chosen as the step 


function discussed above and Pq , is a configuration weighting fupction* 
normalized so /Pg,sin0M6’ = 2^ Since the scattering is independent 


of the laboratory coordinates the differential cross section depends only 
on 3 - (j)’ and not 3 and (f>' so that a change of variables to 
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(j) F (j)* - 3 and 3 3 is warranted. The following integrations may 

2 

be performed: (1) over b which gives a factor R /2; (2) over 3 
which removes 6((J)~3“^) ; (3) and over (p which yields a factor 2 tt 
C ollecting terms Eq. (1.8) becomes 

^ SCXg- VJ d% df, (1.9) 

To perform the integration over 0' the independent and dependent 
variables in <5[Xg-X0i] must be interchanged. This transformation can 
be accomplished by using the relation mentioned above and Eq. (1.9) 
becomes 


cl^ ~ C7-' d I 

I 0-&X. 

which upon substituting Eq. (12) of the text yields 


(I. 10) 


' HIT i\-ZpcC^Y 

This is the "working formula", Eq. 18 of the Sec. III.B. 

The differential cross section including out-of^plane contributions 
may be obtained from (1. 11) by integrating over \p to give 


( 1 . 12 ) 
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APPENDIX II. ROTATIONAL ENERGY CHANGE AS A FUNCTION OF SCATTERIfIG ANGLE 


In this appendix Eq. (2) of Sec. II is solved to find the final 
rotational angular momentum j’ . Once the scattering angle and reaction 
configuration are specified, Eq. (2) may be solved for j' . Note that 
Eq. (2) is a vector equation, implying three independent relationships 
and there are three unknowns corresponding to the components of j' . 

This appendix' is restricted to finding only the magnitude of j’ as a 


function of x • ' ' 

The computation begins by defining the' followihg vectors for a given 
collision configuration; R , from A to the c.m. of BC, r from C to 
the c.m. AB; and r„_ from B to C, along the bond axis. The orbital 
angular momentum vectors are written 

0 ^ L' ^ yU.' 

where y = ACB' + C) (A + B + ,'y' =' (A +‘b)C(A + B + and the 


signs are different because of the definition of § ~ 

gt = w - w- . Specification of r and R requires knowledge of 

the impact parameter. The vector R, is defined in terms the angle 
Y = sin ^(b/R) by 


r R. ( 0 , 


(II. 2) 


and r is found to be 


r ^ - A C^rC) R - f (II. 3) 
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The explicit forms for g and g’ are required for computatipn and are 
given by 


f = fCo>o, o 


(II. 4) 


and 

sXvXs^'Vt ^ ^ (II. 5) 

where ^ is the azimuthal angle ;^or the scattered products. Note for 
this model it is easily shown that = cj) . Since only magnitudes are 
required, Eq. (2) may be squared and written 


= L^ - ^ I 


U' l' + L 


/ 2 . 


(II. 6) 


The cross term j • (L L') has been dropped in Eq. (II. 6) because, 
while it is non^zero, L and L' do not depend on j' and the crqss 
term will vanish when the average is taken over all angles j . Expand- 
ing (II. 6), using the definitions above, one obtains 

Svvv-lS' SiAr^'X. dsvoip 4- QjCr^X (Li“<3 



4 p S tuX dM (^0 - X)+ ^ ^ 




(IIr7) 
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where 

p = -A^A+S>5'R. P ■; -6CA+‘S.+<^^CcA-<B)C6'+c:)3'Vec, 


Equation (II. 7) may now be averaged over all b, ({) , ^nd ijj to obtain 
the angular momentum difference: 

rfi If 


(II. 8) 


The variable 0 is ap impli- 
cit function of ^ is therefore not averaged. When the integrations 

in (II. 8) are performed the following result is obtained; 



where f = C(B + c) This equation is used to compute the results 


discussed in Sec. V.2. 



APPENDIX III. LOADED-^HARD-SPHERE MODEL 


The first model derived as an improvement over the simple hard- 

19 

sphere model for reactive collisions was a model based on the concept 
of loaded-hard-spheres (LHS) . A diagram displaying a reaction configura- 
tion for reactants and defining certain vector quantities is displayed 
in Fig. III.l. 



Fig. III.l Reaction configuration for A + BC. The symbols A, 

B, and C represent the atomic species undergoing 
reaction A + BC AB + C . 

In this figure 6 (6*15) the vector from the c.m. of BC(AB) 

to C(A) and S is a vector between the centers of the two spheres. 
Initially the BC-sphere is loaded; after reaction new spheres may be 
envisioned having the same centers and the same sum of radii but with 
appropriate radii so that AB is now loaded and the load is within the AB 
sphere. The load may appear out of the plane of the page and its orienta-^ 
tion is defined by (B^cj)) as defined in Fig. 1 in the main text. 
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The conservation equations Eqs. (l)-(3) apply and the initial velocity 
g is taken to lie in the positive z- direction. The model definition 
is completed by defining values for three of the ten unknowns in the 
conservation equations. We assume that the momentum transfer occurs as 
it does for hard spheres, along the line of centers, i.e., the impulse 
direction is 



(III.l) 


For the third variable we assume that the final translational energy 
E’ is known or that g i is specified. 

The differences between the LHS model and the model developed in 
the main text (MB model) are that the atom C is at the center of the 
BC-sphere and the direction of the impulse is defined differently. The 
LHS includes an attempt to couple the rotational and orbital angular 
momentum j and L . When a loaded sphere with non-zero | jl travels 
through space it will wobble because it rotates about its own c.m. 

While coupling these variables is an admirable intention it complicates 
the solution of the equations and is probably not worthwhile. In con- 
trast to the MB model the LHS model requires the specification of impact 
parameter , b , to determine the scattering angle X . As can be seen 
from Fig. III.l, because the c.m. of BC does not lie at the center of 

the sphere, the range of impact parameters is a function of , 

'^BC 

This has at least two unpleasant consequences: (1) the limits of 

integration for averaging over configuration are messy, and (2) for 
those configurations in which the BC c.m^ lies above the yz-« plane, the 
function X ” X(b) is double-valued and a rainbow is produced. 


/ 
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The derivation of the LHS equations follows closely the MB deriva- 
tion. Equation (10) is valid here; Eq. (11) is written 

i. ^ (III. 2) 

where (I)z is the z-component of I . The quantity (I)z is identical 
in the two models 

Up to this point the two models agree except that (I)^ cannot be less 
than zero for the LHS. From this point the models diverge, the remain-- 
ing components of I are different and the differential cross section 
is constructed differently. 

Equations (10), (III.l), and (III. 2) and the definition for I 
may be used to find the scattering angle as a function of b depending 
on (6,4)) 


OC 


(KVOX 4 ^ \ ~ 


However, this is not the desired quantity to determine the differential 
cross section, rather X must be the independent variable and b must 
be known as a function of x 


t V) CX ' ©. 


(III. 5) 


From Fig. (III.l) and Eq. (III. 3) one obtains the following equation 
A 


5 '- 




y (III. 6) 
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which gives a quadratic equation in b with the solution 

i. - (s4 i _ (£,,); r ^ 


(III. 7/) 


where 




The positive root is always valid; the negative root is also valid 
(i.e., there exist two b’s for the same x ) when X ^ Xq ^ ^ 30*^, 

The quantity Xq defined by setting b = 0 in Eq, (III. 7) and 

solving for X , 


<^.etv<,Xo= 


(III. 9) 


2 

where D = (6__sin0/S) . Also when c|) > 90®, a rainbow is produced 

BC 


at Xj. determined by 


oce^Xr' F- [Cf-()Cf-‘X'0] 


(III. 10) 


Where F 


2 

= (6^_ sin0 sincf)/S) . The quantity x only defined for 

BC ^ 


(j) > 90°; for <() _< 90° , Xj. “ Xq • Because of the LHS condition 
(I)z >, 0 , Eq. (V.3) implies that for a > 1 there exists a minimum 
scattering angle X^ ^ ® • This condition for a > 0 and a < 0 may 
be written. 






These conditions on the range of x ^^^.y be summarized schematically as 






Z bf<svicW?s _ 'V 


(III. 12) 
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where the range Xq is zero if (j) ^ 90®. 

The differential cross section for the LHS model may be found by 

64 

using the standard formula 


feravi-cWs SvaX 


(III. 13) 


Differentiating Eq. (V.7) with respect to x substituting with 

(III. 13), one obtains 

^ = f^)'{ H(X'lC^)7{CXr~x) ■ 

2 2 

where d cr /d O) is given in Ref. 19. 

KJr 


(III. 14) 


4^ (ocsT 

($ 2.^0 


f (t ^’^ ®<') 7 



(III. 15) 


is defined 

G 


c 


( i - ZsiC^X 



) 


and is the heaviside function which has a value of 1 for positive 
arguments and zero for negative arguments. 

To obtain the final differential cross section, Eq. (JII.13) must 
be averaged over 6 and (j) . This requires dividing 6 and <|) 

into ranges and using limits which are themselves variables of the outer 
integration. When this is done^"^ and the integrations are carried out, 
one obtains the f ar-f rom-obvious result that 
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— — — o (III. 16) 

c\^lO 

Needless to say, this is somewhat of a disappointment I 

The LHS model is considerably more detailed than the simple RP 

19 

model , and thus offers the possibility of more detailed information. 

For example, the final rotational energy as a function of x 
calculated as we have done for the MB model. The calculation is not 
reproduced here (it is more complicated than the MB calculation), how- 
ever, the results are qualitatively the same as those presented in the 
main text for the MB model. A second aspect of the model’s detail is 
the possibility for weighting collisions according to impact parameter 
and/or configuration. When weighting is included the differential cross 
section will depart from the simple RP result. Several weighting 

functions were examined but the findings were discouraging. For all ' 
weighting functions based on configuration, backward scattering is 
still found and the shape of the angular distribution is not radically 
different from the distribution without weighting. Weighting based on 
b yields the same result unless large impact parameters are weighted 
more heavily than small ones, in which case the scattering shifts for- 
ward to side-scattering. It was not possible to obtain forward scatterj[.ng 
This model was therefore abandoned in favor of the MB model. 
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APPENDIX IV. POSSIBLE MODEL MODIFICATION 

The model which has been described is in reality a family of models. 
The original goal was to find a classical model for which the equations 
of motion admit an analytical solution^ This has led to the present 
impulsive model in which the net interaction between the reactants and 
products is compressed into an instantaneous momentum and mass transfer. 
The direction and magnitude of the momentum transfer governs the pre- 
dictions of the model. The Impulse assumption permits the equations 
of motion (or, equivalently, the conservation equations) to be solved 
exactly. A loaded-hard-spheres model was examined first since this 
is the simplest conceivable impulsive model. In this model the momentum 
transfer is constrained to occur along the line-of -centers at collision. 
When this proved unsatisfactory the strict hard-sphere concept was 
abandoned in favor of a reaction sphere determining the point at which 
reaction occurs and the direction of momentum transfer was assumed to 
depend on the collision configuration. The concept and meaning of 
weighting the collision configuration was introduced and discussed. It 
is clear that there is a great deal of leeway here in the development 
of a model. Weighting schemes quite different from the ones suggested 
here can easily be implemented and different criteria for direction of 
the impulse is also possible. In order to solve the equations of motion 
it was assumed that the final translational energy was known. One could 
equivalently specify the final vibrational or rotational energy and obtain 
the final trajectories uniquely. An alternative assumption is to assume 
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that the magnitude of the impulse is known which would then determine 

a final translational energy distribution. A model which assumes that 

the impulse is known and is constant for all collision configurations 

30 

has been developed by Kuntz. 

Many additional refinements may be included which still provide 
an exactly soluble model. Two refinements (which have not been explored 
in depth), directed at improving the final rotational energy distribu- 
tion, are (1) to consider the "spheres" rough so that there is a torque 
perpendicular to the collision axis (this would couple the orbital 
and rotational angular momentum), and (2) introduce the assumption of a 
"breathing sphere," i, e. , define a reaction radius R which is a function 
of the collision configuration. 

For any real system undergoing reaction, the products are not 
produced with one single final translational energy E’ but rather 
with a distribution of final energies which is in general coupled to 
the scattering angle. In order to model a real system, a final energy 
distribution P(E’) or equivalently, if the initial energy is specified, 
a distribution P(ot) , must be guessed (based on say, statistical and 
quantum mechanical consideration) or determined experimentally and then 
folded into the angular distribution. An assumption commonly made in 

this regard is that the angular and energy distribution are uncoupled. 

In the present model this assumption is not involved since the angular 
distribution Eq. (18) depends explicitly on a . The actual dependence 
of Xq ^ is weak, and can easily be modified in one of two ways;: 

either directly by folding in a P(a,y) which depends explicitly on x 
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or indirectly by choosing a weighting function p(6,a) which depends 
on a . Such ad hoc "refinements", however, only serve to weaken the 
value of a model as a predictive tool. They have therefore been dis- 
regarded. 
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APPENDIX V. H + Cs^ 

An interesting case for application of the model is suggested by 
the recent study of the H + CS 2 reaction by Lee, Gordon and Herschbach 
(LGH) . They have approximate results (no velocity analyses were per- 
formed) which indicate the differential reactive cross section is peaked 
in the sideways direction with a half width of - 25®-45®; that it is 
not very sensitive to changes in the initial energy of the reactants; 

and that the most probable final energy is - 1 kcal/mole. Their 
analysis suggested that the H atom has a high probability of attacking 
Cs^ "from the side", causing the anisotropy in the angular distribution. 

Taking a cue from their mechanistic model one should weight most heavily 

those configurations in which Cs^ is perpendicular to the H approach 

vector. (This is a naive abstraction in terms of their rather detailed 

mechanism involving short range HCs - Cs repulsion and long range HCs - Cs 

attraction. However, it seems reasonable to suppose the most probable 

direction of the net impulse is perpendicular to the approach vector if 

this is also the most probable approach.) 

The weighting function chosen to emphasize the sideways attack 

67 

was sin^9 . Because of the high symmetry of this function there is 
no change in scattering when the impulse changes sign so that we are 
searching for a one parameter description (n) of the angular distribu- 
tion (given now, the form of the weighting function and a) . 

As n increases from zero to one the angular distribution 
2 2 

d cr/d tJL) moves from predominently forward scattering to a broad, skewed 
distribution peaking at 60® with a half width of 110®. If the weighting 
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is increased to n = 6 the distribution peaks at 75 “ with a half-width 
of 50°. As n is increased further the distribution remains stationary 
but continues to narrow. To shift the distribution back toward 90°, 
alpha must be Increased. A good fit is obtained for n = 6 . The fact 
that the distribution peaks at 75° is not inconsistent with the data 
reported by LGH. A relatively high value of n is required to give such 
a narrow angular distribution. It may also be significant that the 
wings of the angular distribution change from a convex curvature to con- 
cave, gauss ian-like curvature only for n > 2 . 

The model also predicts, in accord with LGH, that the distribution 
should not significantly change in shape or position as the temperature 
of the reactants is increased. For n = 6, a change of alpha from three 
to twenty shifts the distribution from 65° to 85° and only narrows it 
slightly. Similar results may also be obtained with the model using 
(1 + cos0)’^ weighting with positive impulse for n ~ 0.5. This dis- 
tribution, however, is considerably broadened and for this weighting 
the model predicts a strong temperature dependence of the position and 
breadth of the distribution. Improved experimental information could 
easily distinguish between these two situations. 
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LEGENDS FOR FIGURES 


(a) Momentum diagram illustrating Eq. (6) for the impulse I . 

For the case shown (a < 1), X ^ forward scattering). 

(b) Goordinate diagram, with origin at the center of mass of 

BG . The Initial relative velocity g defines the Z direction. 
R specifies the "reaction radius" and b is the impact para- 
meter. A and the c.m. of BC lie in the XZ-plane . 

Weighting functions R„ , <T„> and <TqP, > . The averages 

U U u b 

are over b , (j) . The functions are arbitrarily scaled to 

unity at 6 = 0 . 

2 2 

. Plots of d a (y)/d 0) for various values of a and for 
K 

R *=» 2 . The weighting function is (1 + cos0)^ where each curve 
is labelled with the choice of sign and the value of n . 

. Plot of configuration angle 6 versus scattering angle x 
for a =0.8. Vector diagrams, similar to Fig. 1, have been 
drawn for small and large values of x * can be seen that 

only configurations with 0 < 90° are possible. 

2 2 

. Plot of d CTj^Cx)/^ ^ for Rb + weighting parameter 

is n = 1 . 

2 2 

. Plot of d cr^(x)/d ^ for the reactions K + Cl^* ^^ 2 ’ ^2 ‘ 

Model parameters: negative impulse, a based on Q = 0, and 
weighting parameter n = 1. 

. Plot of defined in Eq. (21), versus scattering angle x 

for the systems Li, K, Cs + Br 2 * Model parameters are given 
in Table 3. 
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